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In introducing second quantization for fermions, Jordan and Wigner (1927/1928) observed that 
the algebra of a single pair of fermion creation and annihilation operators in quantum mechanics 
is closely related to the algebra of quaternions H. For the first time, here we exploit this fact to 
study nonlinear Bogolyubov-Valatin transformations (canonical transformations for fermions) for a 
single fermionic mode. By means of these transformations, a class of fermionic Hamiltonians in an 
external field is related to the standard Fermi oscillator. 

PACS numbers: 05.30.Fk, 03.65.Fd, 02.10.De 

Keywords: second quantization, Fermi oscillator, isotropic vectors, spinors, coherent states 



Unitary transformations play a prominent role in quan- 
tum mechanics. Like canonical transformations in clas- 
sical mechanics, unitary transformations of quantum dy- 
namical degrees of freedom often simplify the dynamical 
equations, or allow to introduce sensible approximation 
schemes. Such methods have wide-ranging applications, 
from the study of simple systems to many-body prob- 
lems in solid-state or nuclear physics and quantum chem- 
istry, up to the infinite-dimensional systems of quantum 
field theory 0, 0, El- 0] ■ Linear (unitary) canonical trans- 
formations (i.e., transformations preserving the canoni- 
cal anticommutation relations (CAR)) for fermions have 
been introduced by Bogolyubov and Valatin (for two 
fermionic modes) in connection with the study of the 
mechanism of superconductivity d, IE 0, H, E|- These 
(linear) Bogolyubov-Valatin transformations have been 
extended, initially by Bogolyubov and his collaborators 
03, E3> Appendix II, p. 123 (0, p. 116, H, 
p. 679), to involve n fermionic modes [so-called general- 
ized linear Bogolyubov-Valatin transformations, see, e.g., 
[lH - Part III, p. 247 (0, p. 341)]. Such linear canonical 
transformations are important from a physical as well as 
from a mathematical point of view. Mathematically, they 
allow to relate quite arbitrary Hamiltonians quadratic in 
the fermion creation and annihilation operators to collec- 
tions of Fermi oscillators whose mathematics is very well 
understood. From a physical point of view, canonical 
transformations implement the concept of quasiparticles 
in terms of which the physical processes taking place can 
be described and understood in an effective and trans- 
parent manner. To apply the powerful tool of canonical 
transformations to the physically interesting class of non- 
quadratic Hamiltonians, however, requires to go beyond 
linear Bogolyubov-Valatin transformations. Certain as- 
pects of nonlinear Bogolyubov-Valatin transformations 
have received some attention over ti me Il7l|l8l ll9l 12(1 

H3, IlH EH E3, laH laa, 11^,1^1 IHai 

(We disregard here work done within the framework of 



the coupled-cluster method (CCM) -4J which is nonuni- 
tary.). However, a systematic analytic study of general 
(nonlinear) Bogolyubov-Valatin transformations has not 
been undertaken so far. In the present paper, as a first 
step towards this goal we are going to investigate the 
prototypical case of a single fermionic mode. 

Let us consider a pair of fermion creation and anni- 
hilation operators a + , a. Here, we regard the creation 
operator a + as the hermitian conjugate of the annihila- 
tion operator a: a + = at (we will use the latter notation 
throughout). They obey the CAR 

{fit, a} = a* a + a a' = 1, (1) 
(at) 2 = a 2 = 0. (2) 

It is now instructive to consider the following pair of anti- 
hermitian operators. 

a p] = _ a p]t =i(a + at) (3) 
a p] = _ a p]t =a-at (4) 

These two operators obey the equation (p,q = 1,2) 

{aW,a [9] } - -25 pq . (5) 

Consequently, they generate the (real) Clifford algebra 
C(0, 2) which is isomorphic to the algebra of quaternions 
H (cf., e.g., H3, Chap. 15, p. 123, Hi, Chap. 16, p. 205). 
We can define the three quaternionic units i, j, k by the 
equations 

i = a [1] = i(a + a*) , (6) 

J = a [2] = a - ot , (7) 
k = a® = aW &® = i (at a - fifit) . (8) 

Quite generally, these definitions entail that any pair of 
fermionic creation and annihilation operators a', a in- 
duces a (bi-) quaternionic structure into any consideration 
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and model they are a part of. And in turn, any quater- 
nionic structure can be interpreted in terms of fermionic 
creation and annihilation operators. The link between 
the algebra of quaternions H and fermion creation and 
annihilation operators has been mentioned for the first 
time by Jordan and Wiener in introducing second quan- 
tization for fermions |39j, p. 474, [U, p. 635 (H3, p. 45, 
|42j |. p. 113). However, it seems to not have found its way 
into the work of later authors (The only further mention 
of this fact in the literature we have been able to find is 
in ref. |43j.). 

Let us now start by writing down an Ansatz for the 
most general Bogolyubov- Valatin transformation for a 
single fermionic mode. In view of the eqs. JU, (0 
the new pair of fermion annihilation and creation op- 
erators b, ¥ reads (here we assume the coefficients to 
be complex numbers: \^ k '^ 6 C, k, I = 0,1; {A} = 
{A^A^U^A^ 1 )}) 

"b = B({A};fi) 

= A( 0;0 ) + A< 0;1 > a + A^°) at + a^a, (9) 

P = B({A};fi)t 



= A(°;°) + A^ 1 ) at + A( 1 -°) a + AC 1 ; 1 ) fit 6. (10) 
From eq. Q applied to b, w follows: 

2 | A (0;0)|2 + | A (1;0)|2 + | A (0;1)|2 = j (n) 

and from eq. J5J follows: 

(a^) 2 + A^A^ 1 ) = (12) 

while both equations 0J, J5J) also yield 

2A(°;°) + A< 1;1 > = (13) 

([HI, Sect. 2.6, p. 32, eq. (2.91)). Using eq. (O, eq. (TTJ 

can be transformed to read 



| A (1;0)| + | A (0;1)| 



1 



(14) 



(Take absolute values on both sides of the modified 
eq. 112): (A(°;°)) 2 = -A^A^ 1 ), eliminate |A^°)| 2 
from eq. I|ll|) and take the square root.). For compar- 
ison, let us have a look at the class of generalized lin- 
ear Bogolyubov- Valatin transformations (for one mode!): 
A (o ; o) = A (i ; i) = Q Then; eq requires that 

A (i;0) A (0;i) _ o. This condition allows two solutions: 



A d;0) = 0) | A (0;l)| = ^ 
A (0;l) = Q) | A (l;0)| = L 



(15) 
(16) 



It has been found that generalized linear Bogolyubov- 
Valatin transformations (for n modes) are equivalent to 
the group of 0(2n,R) transformations which is in ac- 
cord (for n — 1) with the eqs. (|15f) . (|16fl (This group 



is reduced to SO(2n,IV) if one only allows transforma- 
tions continuously connected to the identity map - then 
in our case only eq. (|T3J) applies: jil IH li a ItR lij 
0113. & Sect - 3 -2, P- 16, Q, Sect. 2.2, p. 38, 
Sect. 9.1, p. Ill, §9.1, p. 127; if one does not as- 
sume that fi and fi + hermitian conjugates of each other 
the correspondi ng groups are 0(2n, C) and S O(2n, C ), 
respectively 0,1110, Q> Sect - 2.1, p. 34, 0,0.). 

The Bogolyubov- Valatin transformation can be in- 
verted. We can write: 



a = B(M;6) 

= + u iO;i) b + u (i;0) g t + V (W) frb. (17) 

Inserting eq. © into eq. i|17|) one obtains a system of 
linear equations in {v} whose (unique) solution reads: 



„(°i°) = A(°;°)A( 1; °) - \<.W)\(W) (18) 

(19) 
(20) 
(21) 



„tti°) = a( 1; °) 

= -2z/ 0;0 > 



One can convince oneself by explicit calculation that the 
{v} given by eqs. (|18|1 - H21|) obey the analogues of eqs. 
(|11|1 . i|12Jl if the {A} obey the latter equations. Further- 
more, the nonlinear Bogolyubov- Valatin transformations 
© form a group Gbv- After the above considerations it 
remains to check that B (-{V}; fi) — B ({/i}; B ({A}; fi)) be- 
longs to Gbv if B ({A}; fi) and B ({//}; fi) belong to Gbv- 
One can explicitly check that the {v} obey the analogues 
of eqs. JTTJ, G2> if the {A}, {^} obey the eqs. JTTJ, ljT2j> . 
or their analogues, respectively. 

To further study the Bogolyubov- Valatin group Gbv 
it turns now out to be useful to consider the linear vector 
space V generated by the operators a, at (V is the space 
of linear operators in Fock space). It is four-dimensional 
and is spanned by the operator basis a T — (1, a, fit, fit a). 
However, taking into account the connection already dis- 
cussed between the operators fi, fit and quaternions it 
turns out to be advantageous to pursue the consideration 
of this linear space in terms of the operator basis (cf. 
eqs. ©-©) a T = (l, fiM, aP\ fi [3] ). The Bogolyubov- 
Valatin transformation @ can be understood as a base 
transformation in the linear space V. We can write 



A({\}) a, 



(22) 



where the 4x4 matrix A ({A}) is a block diagonal matrix 
A = diag(l, A) and A = A ({A}) is the real 3x3 matrix 



/ Re^ ' 1 ' Re^ 1 ' ) Re \ 
Im Im k( 1; °) Im K W 



A ({A}) 



(0;i) 



X K^) X fs^)) X j 



(23) 
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with (taking into account eqs. H11M13|I ) un ^ determi- 
nant (detA = 1) and inverse A({A}) _1 = A({X}) T . 
Here, we have applied the notation: 



K (0;l) = A (0;l) + A (l;0) j 

K (i;0) = i A(o ; i) _ A (l;0) 
K (1;D = A (1;D = _ 2A (0;0) 



(24) 
(25) 

-2K<° i0 >. (26) 



In view of the above considerations the Bogolyubov- 
Valatin group Gbv is equivalent to the group 50(3). 
Given the link between creation and annihilation opera- 
tors and the algebra of quaternions H discussed further 
above this does not come as a big surprise. In accordance 
with eqs. ©-©j the new pair of operators b, w defines 
a transformed system of quaternionic units i', j', k' by 
writing 

i' = &M = i(b + P) 

= Re k (0;1) i + Re k& 0) j + Re k, (27) 
j' = gM = b-P 

= Im k (0;1) i + Im k (1;0) j + Im k, (28) 

k' = m = gm U 2 l (29) 

In terms of the new parameters {k} (eqs. l|24 |l -(|26 |l ) 
the equations l(TT|l . fT^l read 

| K (0;1)|2 + | K (1;0)|2 + ^(1,1) |2 = % ^ 

(f^ + ( K (^°)) 2 + ( K ^ 1 )) 2 = 0. (31) 

Let us now further analyze these equations. Sepa- 
rating them into real and imaginary parts and intro- 
ducing the three-dimensional (complex) vector (e') T = 
(k( 0;1 ), k( 1;0 ), k^' 1 '), these (three real) equations can 
compactly be written as 



(ey e' = 0, |e 



'|2 



e' is an isotropic vector (cf., e.g., |5j|, Sect. 6.3, p. 
113, and for some more detailed and pedagogical 
exposition). In a way, it appears to be an interesting 
feature that within the framework of general (nonlin- 
ear) Bogolyubov-Valatin transformations for a single pair 
of fermion creation and annihilation operators spinors 
make their appearance (via isotropic vectors, cf., e.g., 
[H^|). The properties of these spinors are related to 
canonical (Bogolyubov-Valatin) transformations. Intro- 
ducing two three-dimensional (real) vectors = Re (e'), 
e' 2 = Im (e') (transposed, they agree with the first two 
rows of the matrix (j2HJ)) one can write the eqs. Q32J1 as 



= 1, 



(eife' 2 



= 0. 



(33) 



These equations define the vectors e[, e' 2 as a pair of 
orthonormal vectors which can be supplemented by the 



vector e' 3 = e[ x e' 2 to form an orthonormal vector triple 
in R3. It is worth mentioning here that the vector (e' 3 ) T 
coincides with the third row of the matrix (|23|l . Con- 
sequently, the orthogonality condition(s) for the matrix 
(|23J) are equivalent to the conditions for the Bogolyubov- 
Valatin transformation to be canonical (eqs. H30(l . Ij31|l 
or (|ll(l . (|12Jl ). This is a generalization of an insight ob- 
tained for linear Bogolyubov-Valatin transformation (see 
[JH Oil Oil ) to the general (nonlinear) case. 

The canonical (Bogolyubov-Valatin) transformation 
@ can be implemented by means of an unitary trans- 
formation U ({A}; a): 

b = B({A};a) = U({A};a) aU({A};a) f . (34) 
The analogue of eq. (|3~1|) 

gM = = U({A};a) a [1] U({A};a) f (35) 

has a remarkable interpretation in terms of quaternions 
discussed further above (an analogous comment applies 
to and ai 3 !). Eqs. (|27|l and l|35|l are just concrete 
realizations of the theory of rotations in the language 
of quaternions first elaborated by Cayley and Hamilton 
(cf., e.g., [H, Sect. 12.8, p. 215, eq. (9), 0, Sect. 4.5, p. 
201). Eq. represents a (50(3)) rotation of the vec- 
tor (1,0,0) in the three-dimensional space spanned by 
the quaternionic units i, j, k while eq. (|35() stands for the 
corresponding (5/7(2)) transformation of the quaternion 
i (= fiW) by quaternionic multiplication. The unitary op- 
erator U ({A}; a) can be understood as a unit quaternion 
given by (— it < <f> < it, n±, rii, n?, £ R, n 2 = 1) 



U({A};a) = cos- + sin - (ni i + ni j + n 3 k) (36) 



<f> (m i + n 2 j + n 3 k) /2 



(37) 



The coefficients {A} are given in terms of the parameters 
(32) ) fi\ , ri2, 11,3 by the equations 



A (0;!) — 1 cos _ — m „ am — 
2 6 2 



A (1;0) = (ni + m 2 ) 2 sin 2 

A (l;l) = _ 2A (0;0) 



2' 



(38) 
(39) 



= 2i ^cos — — ins sin — J (n\ + in 2 ) sin — . (40) 

To obtain these relations insert eq. Q36Jl into eq. 134(1 and 
compare the r.h.s. with eq. @. From the representation 
(|37|l one sees immediately that the operators i = aW, 
j = aP', k = fil 3 ! (cf. eqs. ©-©) are generators of the 
group 5/7(2) and they obey the Lie algebra of 50(3), 
SU(2). This has been observed earlier (in a more general 
context) in Ell (also see 0])- Related observations can 
be found in |62|. Appendix A.l, p. 919, [Hjl and [63. p. 
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907, eq. (6.2). One can convince oneself that for linear 
Bogolyubov-Valatin transformations of type l(TB)l eq. i|3T|) 
agrees (sometimes up to some elementary complex phase 
factor) with eq. (7) in jg^, with eq. (5.1) in |ZjJ, with 
eq. (3.6) in HU, with eq. (2.32a), p. 40, Sect. 2.2 in 0, 
with eq. (3.10) in (reduced to the one-mode case; 
incidentally, there is disagreement with [f}?} , P- 205, below 
ofeq. (11)). 

The vacuum state |0) defined by fi|0) = transforms 
under (general, i.e., SO(3)) Bogolyubov-Valatin trans- 
formations according to the law 

|0) {A} = U({A};fi)|0) , 6|0) {A} - 0. (41) 

Associating |0) with a vector in a two-dimensional (com- 
plex) Hilbert space and U ({A}; fi) with a 2 x 2 matrix 
operating in it [cf. [lU, p. 474/475, 0, p. 634 ([Il|, p. 
44, 0] , P- 112)] one sees that this vector transforms as 
a spinor (with a corresponding element of SU(2)) under 
Bogolyubov-Valatin transformations. The state \0){\} is 
a spin (ST/(2)) coherent state [23| (with respect to the 
fi, fit operators, cf., e.g., Ep. Sect. 1.4, p. 25, [5^|. Sect. 
4.3, p. 59, 53], §4-3, p. 72\$g, Sect. III.D.l, p. 884, and 
Sect. VI.A.l, p. 907). However, these fermion coherent 
states are different (cf. the comments in [^|, Sect. 1.5, p. 
55 and in |^, Sect. VI. D, p. 919) from the Grassmann 
(fermion) coherent states (see, e.g., jH^I, Sect. 1.5, p. 48). 

Finally, let us have a look at the standard Fermi oscilla- 
tor given by the Hamiltonian H — fit fi— i. Applying the 
Bogolyubov-Valatin transformation © one can see that 
it is unitarily equivalent to the follo wing class of fermionic 
oscillators in an external field HI El EIH 0, UK, 
Sect. 2.6, p. 29 (below, we have taken into account the 
eqs. ||ni-CElJO: 



H' = 



8tS - ~ = U({A};fi) H U({A};fi)t 



(| A (0;D 



| A (i;0) 



(42) 



As a special case, eq. (|4^|) contains for |A( 0;1 )| = |A( 1; °)| 
also Hamiltonians that are linear in the creation and an- 
nihilation operators (such Hamiltonians have been stud- 
ied in [24, Sect. 4, p. 477). Eq. demonstrates that 
any (Hermitian) Hamiltonian Hq (0 < a £ R, (3 € C) 



a \ a' a 



+ a + /3fiT 



(43) 



can be written and understood in terms of transformed 
creation and annihilation operators , b as 



H = \/a 2 + 4|/^ 



fet b-- 
2 



In view of the above considerations, its dynamical (spec- 
trum generating) algebra is so(3) ~ su(2). 

The present paper paves the way for the study of (non- 
linear) Bogolyubov-Valatin transformations in full gener- 
ality for any finite number of fermionic modes. This is 
done by introducing a methodological framework which 
can be generalized (stepwise) to more than just one mode. 
For several - say n - fermionic modes, the formalism can, 
for example, be expected to allow equivalences of wide 
classes of non-quadratic fermionic Hamiltonians to col- 
lections of n Fermi oscillators to be derived. This will 
be of considerable interest for a wide range of physically 
relevant models. However, beyond its plain methodolog- 
ical value the present study of nonlinear Bogolyubov- 
Valatin transformations for just one fermionic mode pro- 
vides us even with some surprising insight. Note, that the 
nonlinear Bogolyubov-Valatin transformation Q defines 
fermionic operators as a sum of fermi-even and fermi-odd 
terms. This is reminiscent of a supersymmetric transfor- 
mation (cf. in this respect |36lp. In this context, remem- 
ber that linear Bogolyubov-Valatin transformations - in 
contrast to (bosonic) linear Bogolyubov transformations 
- do not allow any linear shifts by complex numbers to 
be performed [see, e.g., |fi£|. Appendix IV, Sect. 1(a), 1. 
ed.: p. 280, 2. ed.: p. 292 T|7b1|. p. 328), also see Q, Sect. 
2.4, p. 40, eq. (2.33b)]. 

The future generalization of the present work to sev- 
eral fermionic modes can alternatively also be understood 
as generalizing the complex coefficients {A} in eq. 10 to 
operator valued functions for which the present analysis 
has to be repeated in an appropriately modified manner. 
Closely related to this direction of future research is to 
consider the coefficients {A} as elements of an appropri- 
ately chosen Grassmann algebra. However, the coeffi- 
cients {A} can not only be imagined to be functions of 
fermionic operators but also to be functions of bosonic 
creation and annihilation operators. Such constructions 
are met, for example, in the study of so-called quantized 
Bogolyubov-Valatin transformations (introduced in [7]]]) 
and, more generally, in the study of boson-fermion inter- 
actions (see, e.g., |2j, Chap. 5, p. 108). 
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